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ABSTRACT 

We extend our previous analysis of the quantum state during and after 0(4)- 
symmetric bubble nucleation to the case including gravitational effects. We find 
that there exists a simple relationship between the case with and without grav- 
itational effects. In a special case of a conformally coupled scalar field which is 
massless except on the bubble wall, the state is found to be conformally equivalent 
to the case without gravity. 



1. Introduction 



Field theoretical quantum tunneling phenomena such as false vacuum decay 
are considered to have played important roles in the dynamics of the universe in 
its early stage. One good example is the so-called extended inflation [1], in which 
the inflationary stage of the universe ends with nucleation of true vacuum bubbles 
and thermalization of the universe by collisions of these nucleated bubbles. 

As another interesting possibility, we have recently proposed a simple one- 
bubble scenario of the inflationary universe by considering the particle creation 
during and after the false vacuum decay [2] . Provided enough entropy is produced, 
it is possible to have our universe inside one nucleated bubble. However, since our 
knowledge of the quantum state after the false vacuum decay is far from sufficient, 
we are unable to argue further for or against this possibility at present. 

Besides the inflationary universe scenario, what happens after the bubble nu- 
cleation is an interesting issue as a fundamental process relating to the quantum 
matter production or the quantification of quantum effects in general in the early 
universe. 

Toward a clear understanding of the issue, various attempts have been made. 
Among them, Rubakov developed a method of non-unitary Bogoliubov transfor- 
mation to treat the particle production during tunneling [3]. Then Vachaspati 
and Vilenkin investigated general features of the quantum state during and after 
nucleation of an (9(4)-symmetric bubble, paying particular attention to the sym- 
metry of the state [4]. Meanwhile, we have developed a method to analyze the 
quantum state during and after field-theoretical quantum tunneling by construct- 
ing a multi-dimensional wave function in a covariant manner [5] (hereafter Paper I). 
Then we have applied it to the (9(4)-symmetric bubble nucleation and investigated 
the properties of the quantum state of fluctuating degrees of freedom in detail by 
constructing an analytically soluble model [6] (hereafter Paper II). 

However, all of these previous analyses were based on several non-trivial as- 
sumptions or simplification of models, the validity of which is not clear. To mention 
one of such, the effect of gravity was neglected in all of them. In this paper, we 
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focus on this point and tackle the problem to incorporate the effect of gravity. More 
precisely, as a first step, we take into account the background spacetime curvature 
induced by the tunneling field solution and investigate its effect on the quantum 
state of fluctuating degrees of freedom. Hence, in particular, the false vacuum is 
de Sitter space. 

This paper is organized as follows. In section 2, we extend our method devel- 
oped in Paper I to the case when the initial state is excited with respect to the 
false vacuum, which is necessary to incorporate the gravitational effect. In section 
3, we show that there is an elegant interrelation between the quantum states after 
tunneling with and without the effect of gravity. As a specific example, we then 
consider a conformally coupled scalar field which interacts with the tunneling field 
on the bubble wall but massless elsewhere, and show that its quantum state is 
conformally equivalent to the case without gravity, i.e., the same result as for the 
Minkowski background obtained in Paper II. However, we also point out a para- 
doxical situation we encounter when evaluating the regularized energy momentum 
tensor. In section 4, we summarize our results. 



2. Multi-dimensional Tunneling Wave Function 

In Paper I, we developed a method to construct the multi-dimensional tunneling 
wave function from the false vacuum ground state in a covariant manner. As we 
shall see in the next section, in order to construct the tunneling wave function in 
curved spacetime, it is necessary to extend the formalism given in Paper I to the 
case when the initial state is in an excited state. Some basic parts of this extension 
have been recently given by Yamamoto [7]. 

We consider a system of D + 1 degrees of freedom whose Lagrangian is given 

by 

L = 1 -G af) (<t>)r^-v(<t>), (21) 

(a,/? = 0,l,. ••,!>; i,j = l,...,D), 

where <p a are the coordinates for the D + l- dimensional space of dynamical variables 
(i.e., superspace) and is the superspace metric. In this section, Greek and Latin 
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indices run from to D and from 1 to D, respectively. For simplicity, we assume 
the potential V{4>) of the form, 

V(<P) = U(X) + ±m* J (X)<j ) ^, (2.2) 

where the tunneling degree of freedom is represented by X = (f)° as a collective 
coordinate and the fluctuating degrees of freedom by <j> 1 , respectively. Further 
we focus on the case when the superspace metric depends only on the tunneling 
degree of freedom, G a = G JX), and has no cross-terms between X and 6 l (i.e., 
G Q . = 0), and assume that the signature of the metric is positive definite. The 
potential U(X) is supposed to have a local minimum at X = X F , which is not the 
absolute minimum, as shown in Fig. 1. We call the point (X, = (X p , 0) the 
false vacuum origin throughout this paper. 

The Hamiltonian operator in the coordinate representation is obtained by re- 
placing the conjugate momenta in the Hamiltonian with the corresponding differ- 
ential operators. In general, there exists the operator ordering ambiguity. Here we 
fix it in such a way that the resulting Hamiltonian takes the form, 

H = -|G^(X)V a V^ + U(X) + \m%(XW4?, (2.3) 

where G a/3 (X) is the inverse matrix of G a ^(X). 

In Paper I, we constructed the quasi- ground- statewave function using the WKB 
approximation. The quasi- ground- statewave function is the lowest eigenstate of the 
Hamiltonian sufficiently localized at the false vacuum. Let us briefly summarize 
the method without rigor. The detailed discussion is given in Paper I. 

(i) First, we impose the WKB ansatz on the wave function, 

* = e 4(ff(»)+w( 1 '+-) ) (2 .4) 

which should solve the time-independent Schrodinger equation, 

HV = E^. (2.5) 

We solve this equation order by order with respect to Ti. The energy eigenvalue E 
is formally divided into two parts, E Q and E v of 0(h ) and O^ 1 ), respectively. 
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The equation in the lowest order of U becomes the Hamilton- Jacobi equation with 
the energy E Q , 

~G a/3 V a W^V p W^ + V(</>) = E Q . (2.6) 

By setting G a/3 V = <p a , Eq.(2.6) gives the Euclidean equation of motion, i.e., 

with respect to the imaginary time r = it. 

(ii) We consider a solution of the Euclidean equation of motion which starts from 
the false vacuum at r = — oo with the zero kinetic energy (i.e., E = E Q := U(X )) 
and arrives at the turning point at r = which is the boundary of the classically 
allowed and forbidden regions. If there are several nontrivial solutions, we choose 
the one which gives the minimum Euclidean action, we call it the dominant escape 
path (DEP), and denote it by 0^(r). In the present case, we have 

(^)(r),0{ o) (r)) =(X(r),0). (2.7) 

(iii) Next, along DEP, we introduce an orthonormal basis, e^(r), lying in the hyper- 
surface S(r) orthogonal to it; G^.(X(r))ea(r)e^(r) = 5 &b , where a runs through 
the range 1, 2, • • • , D. For convenience, we fix the orthonormal basis in such a way 
that e\ := e^r = — oo) diagonalizes := Lu^e^e^, i.e., u &h = u & 5 &h , where 

<4 := m U x )- ( 2 - 8 ) 

Then the required orthonormal basis along DEP is constructed by solving the equa- 
tion, 

®ii + -G ik G h .ei = 0, (2.9) 
which corresponds to a special case of Eq.(2.17) in Paper I. 

(iv) Assuming that W^°\(j)) is known around DEP, we span the hypersurface or- 
thogonal to DEP by the coordinates rj a with respect to the basis e % & and define 



Q , := W<?) 

ab ;ij 



el& h , (2.10) 



where the semicolon represents the (D + l)-dimensional covariant differentiation 
with respect to GLg. We also introduce a matrix K\ which is determined by 
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solving 



dr 2 







K = 


((' 



b • 



(2.11) 



with the boundary condition, 



^W(^Tb> (r--oo), 



(2.12) 



where D/dr and R^ p denote the Lie derivative along DEP and the Riemann tensor 
of the superspace metric G aj3 , respectively. One then finds f2 ab is expressed in terms 
of K , as 

ab 



^ab = ^a C (^ 1 ) cb - 



(2.13) 



(v) With these results in hand, the quasi- ground- statewsve function [5] is found to 
be 



C 



1/4 ! 



OJQT 



2(U(X(r))-E ) J\detK%(r) 



1/4 



(2.14) 



x 



T 

exp (-i J dT '2(U(X(T))-E Q ) + ^TrojT-^Q 



where u Q := U"(X F ), Trcu = X!a=l w a > detcu = []g =1 w a and the normalization 
constant C is given by 



rD 



(2(U(X(r))-E Q )f 4 ro y/4 



C = lim 

r^-oo e2^o r \7Tft/ 



(2.15) 



In the present case, the above result can be re-expressed in terms of the su- 
perspace coordinates (r, (j) 1 ) as follows. First, note that the hypersurface S(r) is 
generally different from the r = const, hypersurface S(r) off the DEP; the latter 
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is warped if G^ ^ 0. Let a point (f, r] a ) on the hypersurface S(f) correspond to 
(r, l ) in the original coordinates. Then from the fact Q(f, ry a ) = Q(t, 4> l ) for any 
scalar function Q and the equation, 



Cl.. = + -G. ., 
y y 2 y 



which follows from the definition of the covariant derivative, where 

d 2 



(2.16) 



^(°) 



i j 4 j ab 



(2.17) 



we find 



(2.18) 



Now, replacing r in Eq.(2.14) with f given by Eq.(2.18), expanding the result 
around r, noting the fact dW^/dr = 2(U(X) — E Q ), and regarding r] a and 0* as 
quantities of 0(h 1//2 ), we find the wave function to the first WKB order to have 
the form, 



*(^) = epo$(x, </>*), 



(2.19) 



where 6 is the lowest WKB part, 



e(x) = 



c 



2(V(X(t))-E q ) 



1/4 



r 

exp (-I | rfr / 2(t/(X(r))-i? ) + ^ r), (2.20) 



and $ is the first WKB correction, 



fdet cu/7r 



1/4 



detX^r)^! 



'(^•- ^H ;/ >v). (2.21) 



Further, from Eqs.(2.13) and (2.16), we obtain 



n.. = VG.i a Xv 



(2.22) 
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Also, from Eq.(2.11), we find K a satisfies the equation, 



G..ki + G..ki = m%Kl, (2.23) 

which is just the classical equation of motion for cj) 1 in the Euclidean time. With 
the help of Eq.(2.23), it is then straightforward to show that to the first WKB order 
(2.21) satisfies 

[V-E[]$ = 0; 

(2.24) 



V 



G-^G 1 / 4 - -G^d.d. + \mM^ 
T 2 1 J 2h 3 



where E[ = hTru. One sees that this is just the Euclidean time Schrodinger 
equation for the fluctuating degrees of freedom. 

Following the procedure taken in [7], we now construct a set of generalized 
creation and annihilation operators, A a and A a whose action on some eigenstate 
of the Hamiltonian produces another eigenstate. In other words, look for operators 
which correspond to the usual creation and annihilation operators at the false 
vacuum origin, i.e., [D, A a ] = w a A a and [D, A a ] = — uJ a A & . To do so, first, we 
make a set of operators a a and a a which commute with the differential operator T>. 
If we assume the forms of a Q and a\ as 




1 ! a <)o' '-'V 2w 



(2.25) 

d 




it is easy to see that a a commute with T> and the necessary condition that a a 
commute with V is that Q a satisfies the same equation as Eq.(2.23) for K a . Then, 
as we have adopted the orthonormal basis e a which diagonalizes at r — > — oo, 
if we set the boundary condition of Q a as 

Qi^e-^ei forr^-oo, (2.26) 

we find [at, oj = 5 &h . Hence the relevant creation and annihilation operators A a 
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and A a are found to be 



e-^a a , 



(2.27) 

An excited state wave function with respect to the fluctuating degrees of freedom 
can be obtained by operating these creation operators, A^, to the quasi- ground- 
statewave function (2.21). 

Given these results, it is convenient to reformulate the method to construct 
the wave function in the following way. We consider a path in the complex plane 
of time as shown in Fig. 2. Along this path, we construct a solution X(t) of 
the classical equation of motion with the initial condition X(— oo) = X p . In the 
segment A of the path, the solution stays at the false vacuum origin X(t) = X p 
which is certainly a solution. This solution can be smoothly connected to the DEP 
solution at a sufficient large negative r (= it). The segment B corresponds to 
the DEP in the forbidden region and C the allowed region. The solution along 
C is obtained by analytically continuing the DEP solution by setting t = —it 
(> 0), which describes the classical motion after tunneling. The lowest WKB order 
wave function is described by this classical solution. Now along this path, we solve 
Eq.(2.23) for K l & and Q l & with the initial conditions at t — > — oo as 



(2.28) 



where u l & and u l * are the positive and negative frequency functions, respectively, in 
the false vacuum. In this way, we obtain the first WKB order wave function before 
and after the tunneling. 

To close this section, perhaps it is worthwhile to mention that the time t (or 
r) discussed here is not the external time in the original Schrodinger equation but 
is a parameter that naturally arises from the characterization of the lowest WKB 
configuration (or the internal time). 
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3. Incorporation of Gravitational Effect 



We consider the system which consists of two scalar fields, i.e., the tunneling 
field a and another field which represents the fluctuating degrees of freedom, 
both coupled to gravity. We consider the situation in which the potential U(a) is 
in the form shown in Fig. 1 and a is initially at the false vacuum, a = a p . The 
interaction between the two fields is assumed to be described by the a-dependent 
mass term of (j); m 2 (a)(j) 2 /2. As in Paper II, we ignore the fluctuations of a and 
the metric g^ v for simplicity. We begin with the Lagrangian of the form, 



C — C qrav + C a + C , , 



(3.1) 



where 



■'grav 



1 



16?rG 



R, 



C 6 = ~ a V^f 



-g^V^V^a + U(a) 



^V^V>+(m 2 (a)+^) 



(3.2) 



and G and R are the gravitational constant and the scalar curvature, respectively. 
Here and in what follows, we use the notation of the 3 + 1 decomposition of the 
spacetime metric: 

-a 2 + (3 s (3 s (3 n 



9 



fim Inm 



(3.3) 



The Hamiltonian is given by 



H = H qrav + H a + H 



(3.4) 



where H grav and H a are the Hamiltonians of the gravitational and tunneling fields, 
respectively, and H, is that of the fluctuating field, 



H, = / d 6 x 



a ^S fc -V^V m ^+(m 2 (a)+^) * 



(3.5) 



with p being the momentum conjugate to <ft. 
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We should note that when gravity comes into play, there exists no external time 
and the Schrodinger equation becomes the Wheeler-DeWitt equation, 



= 0, 



(3.6) 



which is essentially the Hamiltonian constraint for the total system and the su- 
perspace metric has an indefinite signature, presumably with only one timelike 
component. However, discussion of the nature of the Wheeler-DeWitt equation is 



we simply assume that the lowest WKB state is described by a classical solution 
of the (a, g^y) system and ignore problems associated with the Wheeler-DeWitt 
equation. Then, as we shall see below, there arises no conceptual problem with the 
construction of the wave functional for 0. 

3.1 Instanton with gravity 

Let us first construct a non-trivial solution of the Euclidean Einstein-scalar field 
equations, an instanton (or bounce) with gravity, to obtain the lowest WKB order 
picture. In the absence of gravity, it has been shown that the classical solution with 
the minimum action is (9(4)-symmetric [8]. Although it is not proved when gravity 
is present, it seems reasonable that it is also the case in the presence of gravity. 
Hence we assume so. The (9(4)-symmetric instanton with gravity was investigated 
by Coleman and De Luccia [9]. Here we shall not repeat the details but discuss 
only those features of the instanton that will be necessary for our purpose. 

The (9(4)-symmetric instanton takes the form, 



Since we ignore the fluctuations in a and g , we denote the instanton configuration 



out of the scope of this paper. Since we ignore fluctuations in g^ v as well as in a, 




(3.7) 



= 0. 
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simply by a{r]) and g^ijj). Then the Euclidean action becomes 



(3.8) 



where a dot means the derivative with respect to rj. We see that n plays the role 
of a Lagrange multiplier; a consequence of the time reparametrization invariance 
of the system. The variation of S E with respect to a gives 

/A a n\ . dU .„ . 

a+3 * = ( 3 - 9 ) 

\ a n J da 

and that with respect to n gives 

•2 4ttG 2 . 2 2 / 8tvG 2 A 

a —a a —nil — a U (a) I , (3.10) 

which is nothing but the Hamiltonian constraint. Because of this constraint, the 
variation with respect to a does not give an independent equation. 

Let us present the solution of the above equations in the thin-wall case. We 
also assume that the true vacuum energy density U (cr ) is non-negative. We choose 
the gauge 71(77) = a (v) an d indicate the wall position by 77 = 77^. The result is 

a(rj) = 

fa < O, 



H F cosh 77 w 

1 



Hp (cosh(ri - 7] w ) cosh?,^ + sinh(77 - %)y / cosh 2 77^ - (H T /H F ) 2 ^j 
fa > V w ), 



(7(77) = 



a F fa < V)' 

a T fa > V)' 

(3.11) 

where 77^ > always and 

Hl= S -^-U(a F ), H> = ^U(v T ). (3.12) 

A schematic picture of the instanton is shown in Fig. 3. It has the topology of S . 
In the thin-wall approximation, the metric and tunneling field configurations are 
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identical to those of the false vacuum at 77 < rj w , hence at 77 < 0. Here, we restrict 
our attention to the case this holds, but it should be mentioned that this is not 
true in general once the thin-wall approximation breaks down. 

As was pointed out in [10], the coordinate 77 cannot play the role of the 'time' 
parameter r which distinguishes each spatial configuration of the instanton the 
sequence of which connects the false vacuum and the critical bubble configuration 
corresponding to the turning point. A relevant choice of it is obtained by the 
following coordinate transformation, 



cosh 77 



a/ sin 2 t + R 2 cos 2 r 
R 



smr 



cosr 



sin r 



a/ cosh 2 rj — sin 2 r 



a/ sin 2 t + R 2 cos 2 t 
—00 < 77 < 00, < r < 7r), 



R 



smr 
coshry ' 



-7T < T < 7T, < R < I). 



(3.13) 



Then we have 



f dR 2 1 

ds\ = cosh 2 77 a 2 (?7) I (1 - R 2 )dr 2 + - _ r2 + R 2 dtt 2 {2) \ . (3.14) 

How these coordinates span the Euclidean spacetime is schematically shown in 
Fig. 4. 

There are two main differences from the original coordinates. First, these co- 
ordinates reduce to the static de Sitter coordinates when a(rj) = 1 /H p cosh 77. 
Therefore at r < — tt/2, the spatial metric and the tunneling field configuration 
on each r = const, surface are identical to those at the false vacuum origin and 
it is also possible to extend this solution beyond r = — it to t = —00. Second, 
the configuration on the r = surface corresponds to the turning point, where the 
solution can be analytically continued to the Lorentzian region by t = —ir (> 0). 
How these coordinates are analytically continued to the Lorentzian region is shown 
in Fig. 5. Using these coordinates, a full description of the false vacuum decay at 
the lowest WKB order is obtained and our formalism of the multi-dimensional wave 
function can be applied to investigate the next WKB order effects. 
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3.2 Quantum state of the fluctuating field 



Next we turn to the issue of the quantum fluctuations during and after this 
quantum tunneling. As mentioned previously, we neglect the fluctuations of the 
metric and the tunneling field to avoid difficulties. For convenience, we use the 
symbol X to represent the tunneling degree of freedom in the superspace, i.e., 
the spatial configurations of the metric and the tunneling field, = g^(X) and 
a = u(X). We assume X is suitably normalized so that the reduced Hamiltonian 
takes the form, 

H re d --Y^ + U{X)+H <P {X ^ ) - (3 - 15) 

We then investigate the quantum state described by H^(X,<fi) with X now repre- 
senting the DEP parametrized by r, X = X(r). 

To apply our formalism developed in the previous section to the present case, 
we replace the suffices ■ ■ ■ with the spatial coordinates x, y, • • • and a, b, • • • 
with certain eigenvalue indices • • • for a complete set of mode functions, say 




W. \\—K(r)^w k {x,T), 

(3.16) 
V7 (x,t) 

Here one comment is in order. If we wish to deal with the spacetime metric with 
non- vanishing shift vector (3 l , we would need to generalize the formalism in §2 to 
allow the superspace metric to have G Q . ^ 0. However, since there is no shift vector 
in our metric, Eq.(3.14), this generalization is unnecessary at the moment. 

Now as the state before tunneling, the most natural false vacuum state is the 
so-called Euclidean vacuum, which is de Sitter invariant and exhibits the same short 
distance behavior of the field as the Minkowski vacuum [11]. Although the latter 
property is essential to single out the Euclidean vacuum, it can be made explicit 
only when we deal with a specific theory. Hence we focus on the de Sitter invariance 
of the vacuum in this subsection. Since the natural mode functions associated with 
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the form of the metric (3.14) (with a(rj) = (HpCoshr))^ 1 ) does not respect the 
de Sitter invariance, the relation between these mode functions and those for the 
Euclidean vacuum is non-trivial. Specifically, they are related by a Bogoliubov 
transformation. This implies that the Euclidean vacuum is described as some kind 
of an excited state relative to the ground state constructed with respect to the 
static time coordinate (we call the latter the static vacuum for convenience). This 
is the reason why it was necessary to extend our formalism to the case of excited 
states at false vacuum. 

To prepare the Euclidean vacuum and to obtain the Schrodinger wave functional 
relevant for the false vacuum decay, in what follows, we first consider a general static 
spacetime and overview the relation between the usual Heisenberg representation 
of a vacuum state and the corresponding Schrodinger wave functional. Then we 
construct the Euclidean vacuum over the static vacuum and translate the result 
to the Schrodinger picture. Once this is done, it is then straightforward to obtain 
the tunneling wave functional according to the prescription given in the previous 
section. 

Let u^(x,t) be a set of mode functions (not necessarily the positive frequency 
functions with respect to the static time coordinate) and be the corresponding 
annihilation operator; A k \0) = 0, where \0) is the "vacuum" in the Heisenberg 
picture. Then we have 

4>(x, t) = J2 ( u k ( x i t )\ + t) A l) . 

k 

p{x > t] = ^ ^ t)Ak + ^ t)A § ' 

k 



(3.17) 



We assume the mode functions u^(x,t) are orthonormalized with respect to the 
Klein-Gordon inner product, 

(u k ,u p ) := -zjd'x^- (u k ul-u k ul) = 5 kp . (3.18) 

To consider the Schrodinger representation, we introduce time-dependent annihila- 
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tion and creation operators a fc (i) and ot(t), respectively, as 



a 



m (t) = e- iAt A,e iAt , aUt) = e~ iAt AU iAt , (3.19) 



where H is the Hamiltonian operator. The Schrodinger representations of the field 
operators <j> s (x) and P s {x) are given by 

4> s (x) = e- iAt 4>(x,t)e iflt 



( u k ( x ^ t ) a k ( t ) + u *k( x ' t ) a {( t i) > 



k 

l> s {x) t 1111 p{x. I)< '"' ( ' J ' 2()) 



a(a;) fc 



Using these operators, the Schrodinger representation of the vacuum, i.e., = 
e _i-ff *|0) is determined by the condition, 

a k (t)\O(t)) s = 0. (3.21) 

On the other hand, using the orthonormality of the mode functions, aAt) and a\ (t) 
are expressed as 

a k (t) = i f d s x (u* k (x,t)p s (x) - ^u* k (x,t)4> s (x) \ , 
4^ = i J d3x (-\(x,t)p s (x) + ^w fc (a;,t)0 5 (a3)J . 

Then, going over to the coordinate representation by the replacements, 

p s (x) -> -ih-jr^ , s (aj) -> 0(aj) , (3.23) 

we find from Eq.(3.21) that 



|0(f)) s ■ exp [-^r I I d 3 xd 3 y 9.(x. y:t)o(x)q(y) } . (3.2 I) 
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where M is a normalization constant and 

Q{x > y] t] = T^lf ^ ^ t)u * 1{v ' t)j (3-25) 

k 

where uT 1 is defined as 
k 

E ^fc 1 ^* *) = ^ - y)- ( 3 - 26 ) 

Now let us specialize the above to the case of de Sitter space and construct 
the Euclidean vacuum over the static vacuum. Since the Euclidean vacuum is de 
Sitter invariant, we need a set of mode functions which are defined over a complete 
Cauchy surface. However, it cannot be covered by one static chart. Just as in the 
case of describing the Minkowski vacuum in terms of the Rindler mode functions in 
the Minkowski spacetime [12], we therefore need to prepare two static charts. We 
label the quantities associated with these two charts by the indices (1) and (2). The 
two regions are causally disconnected. Furthermore, as clear from the metric (3.14), 
the bubble nucleation takes place only in one of the two regions. For convenience 
we regard the region (1) to be the one in which the bubble nucleation occurs. As 
the time direction is opposite in the two regions, we fix it by identifying the future 
direction with the time direction in the region (1). Thus a complete Cauchy surface 
is given by a hypersurface t = = and the Hamiltonian operator is positive 
in the region (1) and negative in the region (2): 

H, = Hf~Hf\ (3.27) 

where both the operators H^> (i = 1, 2) are positive and have the same form. 
The positive frequency functions in each static chart behave as 

k k 

where and are the static time coordinates in the regions (1) and (2), re- 
spectively. We ort honor malize them by the Klein-Gordon inner product (3.18). 
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(3.28) 



Note that the lapse function a is negative in the region (2). The static vacuum is 
expressed as 

\0) = |0 (1) )®|0 (2) ); A«|O W >=0 (i = 1,2), (3.29) 

where are the annihilation operators associated with the positive frequency 

functions u[ l \ 
k 

The positive frequency functions for the Euclidean vacuum are expressed in 
terms of a Bogoliubov transformation from the static vacuum positive frequency 
functions. Since the regions (1) and (2) are completely symmetric, there exist 
two independent positive frequency functions for the Euclidean vacuum for each 
k, which we denote by and Then in matrix notation with indices k 

suppressed, the Bogoliubov transformation takes the form, 

V 1 ) 



a 















(1) 


a<?>\ 




(2) 


aW ) ' 


= 




= 1, 





' (3M f3 (2) \ (3.30) 



where I denotes the unit matrix. Further, by a unitary transformation, we can 
always make the matrix a diagonal. Hence we may assume 

a kl = a k 6 kk^ a kl = - ( 3 - 31 ) 

Let \0) be the Euclidean vacuum. Then it is characterized by 

A®\O) = (t = l,2), (3.32) 

where are the annihilation operator associated with the positive frequency 
functions u®. The corresponding Bogoliubov transformation for A® is 

Then using the commutation relation between creation and annihilation operators, 
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we find that the vacuum \0) is expressed as 



\0) = Afexp [\ B kk> )A k ]A k< ] ) ® l° (2) )' ( 3 - 34 ) 



2 . 
i,j,k,k 



where A/" is some normalization constant and 



The Euclidean vacuum state is invariant under the action of any generator 

of the de Sitter group. Suppose the indices k of represent the eigenvalues 

k 

associated with spherical mode decomposition, (k,£,m), i.e., each of these modes 
is characterized by an eigenstate of the Hamiltonian H , the angular momentum 
square J 2 , and the ^-component of it J z . Since these operators are generators of 
the de Sitter group, the Euclidean vacuum must be a zero eigenstate of all of them. 
From the fact that all of these operators have the form, Q = — Q^ 2 \ we find 
the matrix B should take the form, 



B =[ b J; b kk> = B k s kk>- ( 3 - 36 ) 

Then from Eqs.(3.31), (3.35) and (3.36), the matrix 8 is found to have the form, 

rvrv rvrv n, r\, r\, 

and 

8* 

B=-^. (3.38) 

k at K ' 

k 

Specifically, the Euclidean vacuum positive frequency functions are given by 



4 2) = vf+« (1) ' 
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and the Euclidean vacuum is given by 



\0) = Afexp B k A k^ A k^ j |0(2) >- (3 - 40) 

At this point, we note that (u[ 2 ^) is proportional to e~ w k l (e* w fe*) for the 

k k 

time coordinate t (= = t^) extended over both regions (1) and (2). This is 
a result of the fact that, although regions (1) and (2) are causally disconnected 
in the Lorentzian regime, they are analytically connected through the Euclidean 
regime as r = it^ and r = ±tt + it^ . It is known that the positive (negative) 
frequency functions for the Euclidean vacuum are characterized by the regularity 
on the upper-half (lower-half) complex t-plane [12]. For later convenience, let us 
consider the negative frequency functions (i.e., the analytic continuations at t — 
and t = —m through the lower-half plane). Then the negative frequency function 
in the region (1), u^* oc e lLU k t( \ is analytically continued to the region (2) as 
til 1 )* oc e 

-u) k TT+iu k t _ Similarly, the negative frequency function in the region (2), 
u^* oc e _la V , is analytically continued to the region (1) as u^* oc e -w fc 7r-Ml, fc* . 
Taking the complex conjugates of them, we then find 

), 

(3.41) 



<2> = a ^> + e-«^>). 



Comparing these with Eq.(3.39) and noting the normalization condition \a, | 2 — 
\P k \ 2 = 1, we obtain 

e i9 k e -w fc ,r+w fc 

a k = , s= , P k = , 9 =, (3-42) 

where 6> fc is a constant phase. Also from Eq.(3.38), we find 

B k = e~"k* . (3.43) 

Now it is easy to interpret the above result to the language of the Schrodinger 
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wave functional. Applying Eq.(3.24) to the present case, we find 

\0(t)) s = Af exp j j d 3 xd 3 y Q(x, y; t)cf>(x)cf>(y)^ , (3.44) 

where 

n(x,y;t) = ^-Yl (^k\x,t)uW-\y,t) + uf{x,t)uf-\y,t)) . (3.45) 
k 

However, as we have seen, this state is a zero eigenstate of the Hamiltonian. Hence 
the wave functional (3.44) should be time-independent. This can be demonstrated 
as follows. Operating e - *^* to Eq.(3.40), we obtain 

|0(t)) s = exp ^ S fc a L 1)t (* (1) ) a i 2)t (* (2) )) \° (1) ( t{1) ))s® \0 (2 \t (2) )) s , (3-46) 

where = = t and ajj)t(t(0) are the time-dependent creation operators as 
defined in Eq.(3.19). Their coordinate representations are given as 

««((«) = J d?x ^ufix^-^-i^-nfix.t^i 

«f (tW) = / d 3 x t »)^ + i^,f ( M)«„)) . 



x) j 

(3.47) 



Since |OW) is the true ground state of the Hamiltonian we have, for = 

= t, 

|0 (1) (*)> s ® |0( 2 )(t)) 5 = (e- iEot \OW)) ® (e^|0( 2 ))) = |0«> ® |0( 2 )>. 

Furthermore, since the Hamiltonian is diagonalized with respect to the static vac- 
uum mode functions, we have 

a^(t)a^(t) = (e-^kUWt) (e^POt) = A^U^ . (3.48) 
Therefore we find 

\6(t)) s = exp {^B k A^A^ |0(D> |0( 2 )). (3.49) 

This form is explicitly time-independent. 
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Now, as we have obtained the asymptotic behavior of the wave functional in 
the false vacuum, we go back to the problem of constructing the tunneling wave 
functional. As clear from Eq.(3.49), the Euclidean vacuum can be considered as a 
superposition of many particle states of the following form, 

(A^) ni ■ ■ ■ (A^) n ™\0^) <g> (A^) ni ■ ■ ■ {A^Y m \0 [2) ). (3.50) 

As there is no causal connection between the regions (1) and (2), the wave functional 
can be constructed independently in each region. First consider the region (1). 
Applying the results obtained in the previous section, the wave functional for the 
region (1) which has the asymptotic form in the false vacuum as 

(A^) ni ■ ■ ■ (A^) n ™ |0W> (3.51) 

fvl rvm 



is given by 



(A^(t^)) m ■ ■ ■ (A^(t^)) nm \0^\t^)) (3.52) 



where \0^\t^)) is the quasi-ground-state wave functional with the parameter 
time aloi 
are given by 



time along the path shown in Fig. 2. The operators A^\t^) and A^(t^) 



4 1 )(t( 1 )) = e-^ (1) a( 1 )(t«) 

f£ K 



= e-^"" / d 3 x (^'(^"^-^"(^'VW) • 
A 0-rt (t (i) )= eVaWtfiW) 

fC fC 



(3.53) 



where «W anc l W W satisfy the field equation, 

fC K 



b^v, - (™V) + ZR)] 4 1} = (* = «> «0» (3-54) 

along the trajectory on the complex t^-plane shown in Fig. 2 with the asymptotic 
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initial condition at < 0, 

w£\x,tM)=uM*(x,tM). (3 ' 55) 

As for the wave functional for the region (2), since the tunneling degree of freedom 

keeps staying at the false vacuum origin, it is ^^-independent. However, we may 

express the wave functional in the same manner as that for region (1) by solving the 

mode functions v[ 2 ^ and w[ 2 ^ along the contour of Fig. 2 on the complex t( 2 )-plane, 
k k 

but with the false vacuum configuration throughout the contour. 

Summing up all the terms again, we obtain the wave functional which describes 
tunneling from the Euclidean vacuum, 

$ [(/)( •);*]= A/"exp (j2B k A^(t)A^ \0M(t)) ® \0^). (3.56) 

Now, as we have discussed when deriving the Bogoliubov coefficients (3.42), 
when we consider the negative frequency functions, the regions (1) and (2) are 
analytically connected through the lower-half complex t-plane while the tunneling 
field is at the false vacuum origin. Hence we expect the tunneling wave functional to 
be obtained by finding the mode functions for the tunneling background which 
correspond to the Euclidean vacuum negative frequency functions u^*. Namely, 
first we set the boundary condition, 

w®(x,t) = u®*(x,t), (3.57) 

at t < 0. Then they solve the field equation (3.54) in the Lorentzian time to 
t = = = 0, and further in the Euclidean time with r = it^ and r = 
— 7T + it^ beyond r = — n to r = — oo. Note that from Eq.(3.41), the solutions 
are automatically consistent with the analytic continuation at both r = and —n. 
Then solving back to r = through the non-trivial 0(4) bubble background only 
in the interval —tc/2 < r < 0, and analytically continue to the Lorentzian time 
with t = t^> = — z(r + 7r) at r = — n to the de Sitter background and t = = —ir 
at t = to the (9(3, 1) bubble background. 
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Then, with thus obtained mode functions w[ l \ we can in fact show that the 

k 

tunneling wave functional satisfies 

afm[<P(-);t] = 0, (3.58) 

where 

fiW(f) = J d*x (hw^^x^^-i^fix^x)^ . (3.59) 

After all, as clear from the above procedure, we do not have to know «)W at r < — n 

k 

but only those in the interval — n < r < 0. Furthermore coincide with «W 
J - - k k 

in the region r < —tt/2. Once we become aware of these facts, we need not stick 

to the construction of themselves, nor to the coordinates of the metric (3.14). 

k 

Instead, any complete set of mode functions which are related to u)W by a unitary 
transformation is relevant and any convenient coordinate system can be chosen to 
solve for them. For this reason, we may drop the superscript (i) for the mode 
functions and denote them simply by w^. Thus the procedure to construct the 
tunneling wave functional in the case with gravity turns out to be very similar to 
the case without gravity. In particular, the resulting quantum state after tunneling 
will be related to the true vacuum state by a Bogoliubov transformation as discussed 
in Papers I and II. 

3.3 CONFORMAL SCALAR MODEL 

Here, as a simple application of our formalism, we consider a conformally cou- 
pled scalar field <j> (i.e., £ = 1/6 in Eq.(3.5)) which is massless except on the 
bubble wall. In Paper II, we have investigated a similar model in the absence of 
the background curvature. We show below that this conformal scalar model gives 
the quantum state which is conformally equivalent to the one without gravity in 
Paper II. 

A convenient choice of the coordinates for the present case is obtained by the 
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coordinate transformation of the metric (3.7) with 77.(77) = a (v) as 

T = —e r ' cosr, 
p = e 71 sin r, 

or equivalently that of the metric (3.14) as 

Vl-R 2 sin r(l + Vl-R 2 cos r) 



(3.60) 



T 



E sin 2 r + .R 2 cos 2 r 

i?(l + Vl - R 2 cost) 
sin 2 r + R 2 cos 2 r 



(3.61) 



Then the metric becomes 

d8% = Sl%(t E ) (dT 2 + dp 2 + p 2 dn 2 ^ , (3.62) 



where 

„ a 2 (ln£„) / 

"1 = -^; t E =jT 2 +p 2 . (3.63) 



As clear from Eq.(3.61), the analytic continuation to the Lorentzian metric by 
t — it (t — t^) and r = — n + it (t = t^) corresponds to that at T E = by 
T E = iT. Hence the Lorentzian version of the metric (3.62) is 

ds 2 = fi 2 (0 (-dT 2 + dp 2 + p 2 dQ 2 {2 ^j , (3.64) 

where 

n 2 = ^^; ^V^T^T?. (3.65) 
At the false vacuum origin, Eq.(3.64) reduces to the de Sitter metric: 

d*Ls = ^desiO (-dT 2 + dp 2 + p 2 dn 2 2) ) ; 

2 (3-66) 

^deS = H F (l-T 2 + p 2 ) ■ 

As we have seen in the previous subsection, the procedure to obtain the tun- 
neling wave functional is to solve for the mode functions in the interval — it < 
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t < of the Euclidean tunneling background and analytically continue them to 
the Lorentzian background at r = —n and r = with the condition that are 
unitarily equivalent to the analytic continuation of the Euclidean vacuum negative 
frequency functions in the region r < —tv/2. Hence we may regard Eq.(3.64) 
to represent the background metric both in the Euclidean and Lorentzian regimes 
by allowing T to take complex values. We show how the time coordinate T spans 
the spacetime in Fig. 6. The advantage of using the coordinates (T, p) is that the 
correspondence to the flat spacetime case becomes transparent. The metric (3.64) 
is conformally equivalent to the flat spacetime metric: g^ v = Vt 2 rj^ v . Further, the 
conformal factor Q is a function of only £, i.e., it is 0(3, 1)- (or 0(4)-) invariant. 
In particular, Q may be regarded as a function of the tunneling field a. 

Because of the conformal coupling of the scalar field, the field equation (3.54) 
for the mode functions can be conformally transformed to that on the flat 
spacetime: 

[iTV^ - {m\a)n\a))] w fk = 0, (3.67) 
where we have regarded Q as a function of a and 



w fk = Qw k . (3.68) 

Thus we can construct the mode functions which satisfy Eq.(3.54) by solving 
Eq.(3.67) in the flat spacetime. 

The remaining task is to impose the correct boundary condition on w ... For 

JK 

this purpose, let us consider the case of pure de Sitter background. If the conformal 
vacuum defined by the positive frequency functions u , (oc e~ lkT ) on the flat space 

JK 

agrees with the Euclidean vacuum, the required initial condition for is trivial; 

Wj, k = u *f K - To see this is indeed the case, we examine the symmetric two-point 

function in the conformal vacuum. It is well-known that the positive frequency 

function u. for this conformal vacuum is given by u. = and as a result 

k ° J k fk 

the two-point function is given by 

GW(x,x') = n-} s {x 2 )D^\x,x')n-} s {x% (3.69) 
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where & deS is given by Eq.(3.66), x 2 = -(x ) 2 + x 2 = -T 2 + p 2 and 



B(1) <^'> : = (3 ' 70) 



is the symmetric two-point function in the Minkowski vacuum. 

If we embed the de Sitter space in the five dimensional Minkowski space, 



ds 2 = -(dx ) 2 + (dx 1 ) 2 + (dx 2 ) 2 + (dx 3 ) 2 + (c£r 4 ) 2 , (3.71) 

we find the coordinates x a (a = 0, 1, 2, 3, 4) on the de Sitter space are expressed in 
terms of x^ = (T, x) as 



It is then easy to show that 

° m ^ = hwh?- (3 - 73) 

This shows is de Sitter invariant. Hence, together with the fact that its short 
distance behavior is the same as that in the Minkowski vacuum, the present con- 
formal vacuum is found to be the Euclidean vacuum. The relevant mode functions 
for the tunneling wave functional are then given by solving Eq.(3.67) along the 
contour shown in Fig. 2 on the complex T-plane and multiplying the result by the 

inverse of the conformal factor; w = Vt^w 

k j ft 

To summarize, in the present case of a conformally coupled scalar field, the 
effects of gravity to the quantum state after tunneling is to solve the mode functions 
for the flat background with the mass term m 2 replaced by m 2 Vt 2 and multiplying 
the resultant mode functions by with VL given in Eq.(3.65). The two-point 
functions for the quantum state after tunneling are also given by the conformal 
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transformation of those obtained for the flat background D(x, x') as 



G(x, x') = tt-\x 2 )D(x, x')n-\x' 2 ). (3.74) 

However, the evaluation of the energy momentum tensor seems to require some 
care. It is known that the regularized vacuum expectation value of the energy mo- 
mentum tensor, (T^ u ), for a conformally coupled field on a conformally flat space- 
time consists of term arising from a trivial conformal transformation of ( T^ v ) 

\ J / reg 

in flat space and the terms representing conformal anomalies [13]: 

(T fIU ) reg = f2~ 2 (rlvj + conformal anomalies. (3.75) 

It is also known that if there arises no further divergence apart from the common 
ones for any spacetime, Eq.(3.75) continues to hold for arbitrary state. However, in 
the present case, as discussed in Paper II, we encounter a new type of divergences 
which may be partly due to the delta-function nature of the mass term in our model 
and also possibly due to the breakdown of the WKB expansion. Unfortunately at 
the moment, we are unable to clarify if these new divergences would give rise to 
terms which are conformally non-trivial. If not, since the finite terms of 
we have found in flat space in Paper II will be absent in the case of conformal 
coupling, except for the term which diverges on the light cone, we would find 
only the conformal anomaly terms in the present model, provided we impose the 
regularity on the light cone. Furthermore, if the true vacuum has no vacuum energy, 
it reduces to the flat space and we would find no finite term at all, which sounds 
rather paradoxical. In Paper II, we have argued that the regularity on the light 
cone is necessary to keep the validity of the WKB expansion. Hence if we allow 
the presence of the term which diverges on the light cone, it will be necessary to 
seriously consider the possible breakdown of the WKB expansion. The resolution 
of this issue is left for future study. 
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4. Conclusions 



We have considered an extension of our previous analysis in Paper II of the 
quantum state after 0(4)-symmetric bubble nucleation in flat space to the case 
with the gravitational effect. In order to do so, we have first extended the formalism 
developed in Paper I to the case of multi-dimensional tunneling from an excited 
state at the false vacuum origin. Then using the result of extension, we have 
developed a method to obtain the tunneling wave functional from the false vacuum 
to the true vacuum through a non-trivial geometry of the background spacetime 
described by the 0(4)-symmetric bubble with gravity. Provided that the 0(4) 
bubble is described by the thin-wall approximation, we have found the procedure 
to construct the tunneling wave functional can be formulated in a quite similar 
manner as in the case of flat spacetime background. 

As an explicit demonstration of our formalism, we have considered a simple 
conformal scalar model which is massless except on the bubble wall to represent 
the fluctuations around the 0(4) bubble. We have then found the the resulting 
quantum state is conformally equivalent to that in the absence of gravity, i.e., it 
is described by a Bogoliubov transformation of the true vacuum state (a squeezed 
state). However, we have argued that the evaluation of the regularized expectation 
value of the energy momentum tensor for this quantum state may be highly non- 
trivial, apart from the conventional conformal anomalies. We have also pointed 
out the paradoxical situation that the regularized energy momentum tensor might 
vanish due to the conformal coupling nature of our model if the true vacuum has 
no vacuum energy density. 

At the moment, we are unfortunately unable to judge whether these issues are 
particularities of our over-simplified model or intrinsic difficulties associated with 
field theoretical tunneling phenomena. Further research on the present subject is 
apparently required. 



- 29 - 



ACKNOWLEDGEMENTS 

This work was supported by Monbusho Grant-in-Aid for Scientific Research 
Nos. 2010 and 05640342, and the Sumitomo Foundation. 

REFERENCES 

1. D. La and P. J. Steinhardt, Phys. Rev. Lett. 62 (1989) 376; 

P. J. Steinhardt and F. S. Accetta, Phys. Rev. Lett. 64 (1990) 2740; 

F. C. Adams and K. Freese, Phys. Rev. D43 (1991) 353; 

R. Holman, E. W. Kolb, S. L. Vadas and Y. Wang, Phys. Lett. B269 (1991) 
252; 

A. R. Liddle and D. Wands, Phys. Rev. D45 (1992) 2665. 

2. M. Sasaki, T. Tanaka, K. Yamamoto and J. Yokoyama, Phys. Lett. 317B 
(1993) 510. 

3. V.A. Rubakov, Nucl. Phys. B245 (1984) 481. 

4. T. Vachaspati and A. Vilenkin, Phys. Rev. D43 (1991) 3846. 

5. T. Tanaka, M. Sasaki and K. Yamamoto, Phys. Rev. D49 (1994) 1039, 
[Paper I]. 

6. M. Sasaki, T. Tanaka, K. Yamamoto and J. Yokoyama, Prog. Theor. Phys. 
90 (1993) 1019, [Paper II]. 

7. K. Yamamoto, Prog. Theor. Phys. 91 (1994) 437. 

8. S. Coleman, V. Glaser and A. Martin, Commun. Math. Phys. 58 (1978) 211; 

S. Coleman, in The Whys of Subnuclear Physics, ed. A. Zichichi (Plenum, 
New York, 1979), p.805. 

9. S. Coleman and F. De Luccia, Phys. Rev. D21 (1980) 3305. 

10. T. Tanaka and M. Sasaki, Prog. Theor. Phys. 88 (1992) 503. 

11. B. Allen, Phys. Rev. D32 (1985) 3136. 

- 30 - 



12. W. G. Unruh Phys. Rev. D14 (1976) 870; 

N.D. Birrell and P.C.W. Davies, Quantum fields in curved space, (Cambridge 
University Press, Cambridge, 1982). 

13. P. Candelas and J.S. Dowker, Phys. Rev. D19 (1979) 2902. 



Figure Captions 

Fig. 1: The potential form of the tunneling field, where a p and a T represent the 
values of a in the false and true vacua, respectively. 

Fig. 2: A path on the complex plane of time, which represents a tunneling process. 

The segments A, B and C correspond to the motion staying in the false 
vacuum, an instanton and the motion after nucleation, respectively. 

Fig. 3: A schematic picture of the Coleman DeLuccia instanton solution in the thin 
wall limit. Two dimensions are suppressed 

Fig. 4: Foliation of the Coleman DeLuccia instanton solution by r = const, hyper- 
surfaces. 

Fig. 5: The analytic continuation to the Lorentzian region of the Coleman DeLuccia 
instanton solution. The lower and upper halves are Euclidean and Lorentzian 
regions, respectively. 

Fig. 6: The same as Fig. 5, but with foliation by the T = constant hypersurfaces. 
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